TeV 7 rays from photodisintegration and daughter deexcitation of cosmic-ray nuclei 
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It is commonly assumed that high-energy 7-rays are made via either purely electromagnetic processes 
or the hadronic process of pion production, followed by decay. We investigate astrophysical contexts 
where a third process (A*) would dominate, namely the photo-disintegration of highly boosted 
nuclei followed by daughter de-excitation. Starbust regions such as Cygnus OB2 appear to be 
promising sites for TeV 7-ray emission via this mechanism. A unique feature of the A* process 
is a sharp flattening of the energy spectrum below ~ 10TeV/(T/eV) for 7-ray emission from a 
thermal region of temperature T. We also check that a diffuse 7-ray component resulting from the 
interaction of a possible extreme-energy cosmic-ray nuclei with background radiation is well below 
the observed EGRET data. The A* mechanism described herein offers an important contribution 
to 7-ray astronomy in the era of intense observational activity. 



Introduction. In the field of TeV 7-ray-astronomy, 
new instruments are discovering new sources at a rapid 
rate, both within our Galaxy and outside the Galaxy 
Not surprisingly, one of the brightest TeV 7-ray sources 
is the one discovered first, the Grab pulsar wind neb- 
ula. The integral 7-ray flux obtained from the Crab by 
the Whipple GoUaboration is now the standard TeV flux 
unit: fcrab(e!;'^^ > 0.35 TeV) = lO-^VcrnVs 0. The 
spectral index of the Grab's integrated flux is measured to 
be —1.5, so Fcrab falls by a factor of 30 for each decade of 
increase in min(ejj;^^). Newly commissioned atmospheric 
Gerenkov telescopes (GANGAROO, HESS, MAGIG and 
VERITAS) will reach a sensitivity 100 times below Fcrab- 

Two well-known mechanisms for generating TeV 7- 
rays in astrophysical sources 3] are the purely electro- 
magnetic (EM) synchrotron emission and inverse Gomp- 
ton scattering, and the hadronic (PION) one in which 
7-rays originate from tt" production and decay. There 
has been considerable debate over which of these two 
mechanisms is dominant. The mechanism for the PION 
mode may be either pp or p^ collisions. In this Let- 
ter, we highlight a third dynamic which leads to TeV 7- 
rays: photo-disintegration of high-energy nuclei, followed 
by immediate photo-emission from the excited daughter 
nuclei. For brevity, we label the photo-nuclear process 
A + 7 A'* + X, followed by A'* — > A' + 7-ray as 
"A*". It is likely that each of the three mechanisms 
is operative in some environments. As we show below, 
the A* -process is likely operative in the interesting com- 
plex environments which have accelerated nuclei and hot 
starlight. Fortunately, future measurements can read- 
ily distinguish among the three competing mechanisms 
per source. In particular, the nuclei of the A* process 
act in analogy to a relativistically moving mirror Q, to 
"double-boost" eV starlight to TeV energies for a nuclear 
boost factor F^ = E'^^^/to^ > 10*^. 

The A*-Process. Nuclear photo-disintegration in the 



astrophysical context was first calculated by Stecker in a 
seminal paper almost forty years ago. To our knowl- 
edge, the possibly important role of nuclear de-excitation 
in the astrophysical context was first appreciated and 
proposed by Moskalenko and collaborators more than 
a decade ago. Since then, the A*-process has been over- 
looked by the 7-ray community. Now that data is becom- 
ing available which can validate or invalidate the process, 
it is timely to revive and further develop the A*-process. 
We do so, by providing calculational details, and by iden- 
tifying the astrophysical context where this process might 
dominate over the EM and PION modes for production 
of 7-rays. A detailed discussion on the main features of 
the EM and PION mechanisms is presented in a longer 
accompanying paper Q . 

By far the largest contribution to the photo-excitation 
cross-section comes from the Giant Dipole Resonance 
(GDR) at e™^ ~ 10 MeV - 30 MeV in the nuclear rest 
frame fsl . The ambient photon energy required to excite 

/Ta. The GDR decays 
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the GDR is therefore e 
by the statistical emission of a single nucleon, leaving an 
excited daughter nucleus {A — 1)*. The probability for 
emission of two (or more) nucleons is smaller by an or- 
der of magnitude. The excited daughter nuclei typically 
de-excite by emitting one or more photons of energies 
1 — 5 MeV in the nuclear rest frame. The lab- 
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frame energy of the 7-ray is then e. 
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As just outlined, the boost in the nuclei energy from 
rest to E^^^ = TaAitin plays two roles. It pro- 
motes the thermal energies of the ambient photons to 
the tens of MeV 7-rays capable of exciting the GDR, 
and it promotes the de-excitation photons from a few 
MeV to much higher energies, potentially detectable in 
7-ray telescopes. Eliminating Ta above leads to the 
relation eeh^'^ - eS^^e^''". The "~" indicates that 



the 
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and 



^dxn 



energies are not sharp, but rather 



distributed over their resonant shapes ( "Lorentzian" or 
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"Breit- Wigner" ) . Each of the two spectra have a width 
to mass ratio less than unity. It is sufficient for our 
purposes to treat each spectrum in the narrow-width- 
approximation (NWA). Well-defined values for e^°^ and 
gdxn ^jjgjj result. We may summarize up to this point 
by saying that the A* process produces ^ -ray's with en- 
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exists an accelerated nuclear flux with boost Ta 
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7 X 10^ (T/eV) or equivalent energy 



LAB 



7PeV/(T/eV)pernucleon [6|,l9l. 

Here and throughout we assume a Bose-Einstein (BE) 
distribution with temperature T for the ambient photons; 
in the mean, (e) ~ 3 T. As an example, 7-ray's with 
energy e}^^^ ~ 20TeV are generated in this A* -process 
if an ambient photon temperature of an eV and a boost 
factor of 7 X 10^ are present. Thus, we have arrived at 
an astrophysical environment where the A*-process may 
dominate production of TeV 7-ray's: the region must 
contain far-UV photons (commonly defined as 1-20 eV) 
from the Lyman-a emission of young, massive, hot stars 
such as O and B stars which have surface temperatures 
of - 40, OOOK (3.4 eV and (e) ~ 10 eV) and 18,000K 
(1.5 eV and (e) ~ 5 eV), respectively; and the region 
must contain shocks, giant winds, or other mechanisms 
which accelerate nuclei to energies in excess of a PeV per 
nucleon. Violent starburst regions, such as the one in 
the direction of Cygnus, are splendid examples of regions 
which contain OB stars, shocks and giant stellar winds. 
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fixed at their central values, the 



7-ray spectrum dn[e}^^^) / de}^^^ is given by a simple Ja- 



cobian times the BE distribution with argument e set to 
i.e.. 
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Three regions in e!j;^^ emerge. For e!!j;^^ < e™ ^ £^ ™/T, 



the 7-ray spectrum is exponentially suppressed [11|. For 
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for the A* process, but it is an order of magnitude larger 
at rriTr for the PION p7 process [l^. This means that 
for fixed ambient temperature, need be an order of 
magnitude larger for the PION p7 process, and the re- 
sulting 7-ray energies, proportional to F^, are two orders 
of magnitude larger. The EM and PION pp processes 
contrast with the A* -process in that there is either no 
threshold (EM) or very small threshold 0(2to^) in the 
lab (PION pp), and the resulting 7-ray spectrum rises 
monotonically with decreasing e^"^^ . 

The A*-Rate. We now derive the rate for 7-ray pro- 
duction in the yl*-process. The cross-section is domi- 
nated by the GDR dipole form [l^l, which in the NWA 
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where r^^^ and cr*^^^ are the GDR width and cross- 
section at maximum. Fitted numerical formulas are 
^GDR ^ ;l 45^ X 10-27cni2, pC^DR = 8 MeV, and 
fGDR ^ 42.65^-"-2i for ^ > 4 and 0.925^2-433 

for < 4 [l7|. The 8 MeV width implies a very short 
de-excitation distance of ^ 25 Ta fm. 

The general formula for the inverse photo- 
disintegration mean-free-path (mfp) [l8| for a highly 
relativistic nucleus with energy Fa/to^ propagating 
through an isotropic photon background with energy e 
and spectrum dn{e)/de is [1, 0] 
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For a nucleus passing through a region of thermal pho- 
tons, integration over the BE photon distribution gives 
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gdxn^rp^ there is power-law suppression. where we have defined a dimensionless scaling variable 



The 7-ray spectrum peaks near e: 
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/T. In 
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/2 Ta T. From the prefactor, we learn that the 



particular, one notes that the suppression of the high- 
energy Wien end of the thermal photon spectrum has 
led to a similar suppression of the photon spectrum below 
e^^^ ^ 20TeV/(T/eV) . This lower-energy suppression 
presents a robust prediction of the A* -process. Moreover, 
it contrasts greatly with the PION and EM processes, 
and so provides a unique signature. The A*process 
predicts "orphan" sources, with suppression of associ- 
ated GeV 7-ray's or MeV X-rays (although the photo- 
dissociated neutrons may /?-decay to neutrinos [l3|)- Ob- 
servationally, orphan 7-ray sources have been identi- 
fied [13, [ij], and some orphans are known to be near 
OB starburst regions [14 1. 

Comparing to the PION p7 process, the A* process has 
a much lower ambient photon energy threshold; in the nu- 
clear rest frame, the photon threshold is ~ 10 MeV 



peak of {Xa^)~^ scales in A as cr^°^/e, 
the value of F^i at the peak scales as e. 
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The scaling function f{w) 
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shown in Fig. [T] Approximations to the | In | term yield 
e^™ for w > 2, and | lnw| for w ^ 1. Thus, the exponen- 
tial suppression of the process appears again for large for 
w > 2, i.e., for e^AB < ^gdr gdxrY4 ^^^^ ^ 

region presents a mfp that scales as | lnw|. The peak 
region provides the smallest inverse mfp, and so this re- 
gion dominates the A*-process. In the peak region, w is 
of order one, which implies that Ta T ^ e^^^. When this 
latter relation between the nuclei boost and the ambient 
photon temperature is met, then the photo-disintegration 
rate is optimized. 

The area under the peak region in Fig. [T] is of order 
one, which leads to a simple and reasonable estimate of 
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FIG. 1: The scaling functions f(w) and Xf^F{xt), propor- 
tional to the photo-dissociation rate in Eq. (4) and the 7-ray 
dF/ de spectrum in Eq. (8) , respectively. 



the inverse mfp given in Eq. ^ , for any nucleus boosted 



near Ta 
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/T. The useful and eminently sensible 



estimate is (A^ ) 
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the typical value , 
bers, this estimate yields 



here, we have input 
1/3. Putting in num- 
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5 X 10^^ cm 



3AU 
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for a nucleus with energy in the peak regime around 



10yl/(T/eV) PeV. Here we have used 
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Ea 

2C(3)T3/7r2 ~ T3/4. 

An important question is how many photo- 
disintegration steps in the nuclear chain A ^ [A — 1)* 
[A — 2)* etc. may be expected. Each step will produce 
an excited daughter which may de-excite via emission 
of 7-ray's, each having a typical lab-frame energy 
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^GDR ^dxn 



e^J™/3T. Clearly, the number of steps 
depends on the mfp of each excited daughter nucleus, 
and on the diffusion time of the nuclei in the thermal 
region. The probability for n-cascades, producing n 
excited daughter nuclei and n single nucleons, may be 
written in symmetric ways as 
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(6) 

where = 0, Xn+i equals the diffusion length D of nu- 
clei in the A*-region, the ordered xi < X2 ■ ■ ■ Xn denote 
the spatial positions of successive photo-disintegrations 
to (^-1)*, ... {A-n)*,6j = Aj^-Aj+i, and An+i =00. 
The exponentials in ^ are probabilities that the various 
daughters not interact (i.e., survive) from the point of 
creation to the point of photo-disintegration. A simple 



result obtains when the mfp's are long on the diffusion 
scale D of the ^*-region. In this case the nested inte- 
grals collapse to P„ w ("0^^ ^^=1 E>/Xj- This is just 
the product of the independent probabilities for each ex- 
cited daughter to be produced, times the factor 1/n! that 
divides out all but the one correct time ordering of the 
n photo-disintegrations. One sees that, since D/\j <^ 1 
by assumption, disintegration to more than the first ex- 
cited daughter is unlikely. Such is the case if photo- 
disintegration of nuclei occurs in the Cygnus OB star- 
burst region, as we now demonstrate. 

Relevance to Starburst Regions. The aforemen- 
tioned possibility of generating Galactic TeV 7-ray's from 
accelerated nuclei scattering on starlight in starburst re- 
gions is of considerable astrophysical interest. The energy 
of the photons (~ ZTi, in the mean) is maintained as the 
photons disperse from the stars, but the photon density 
decreases by the ratio of the surface of stellar emission 
(A'i X 4:11 R^) to the loss surface of the starburst region 
(47ri?gg). Taking the giant star radius to be i?^ ~ 10 i?©, 
the radius of the starburst region to be i?sB ~ 10 pc, and 
the stellar count to be A^* ^ 2600, the photon density is 
diluted by 10^^ According to Eqs. ^ and dH), (A^)"^ 
is proportional to the photon density, and hence to this 
factor. Including this factor in Eq. ((S]), one gets the esti- 
mate Xa ^ (56/A) (1.5 eV/T^)^ x 10^^ cm for the nuclear 
mfp in the starburst region. The hot stars are domi- 
nantly (95%) B type, with T^ ~ 1.5 eV. We note here 
that the diffusion time for a nucleus in the Cygnus OB 
starburst region is calculated to be ~ lO** yr ~ 10^^ cm, 
an order of magnitude shorter than the nuclear mfp. Ac- 
cordingly, one expects only a few percent of the nuclei to 
photo-disintegrate in the Cygnus OB region, with multi- 
ple disintegrations being rather rare. 

Some early statistical-model calculations for the pro- 
duction of 7-rays through the decay of the GDR gave 
a mean photon multiplicity between 0.5 and 2 for the 
different nuclides in the chain reaction (see Ref. Q for 
details). We will simplify the data by assuming that one 
photon is emitted per nuclear de-excitation. Then, the 
rate of photo-emission is just the rate of excited daugh- 
ter production, which in turn is just the rate of photo- 
disintegration. Allowing for the 1/r^ dilution of the 7-ray 
flux from the source region of volume Va* at distance d, 
the observed integral 7-ray flux T-y at Earth with e^"^^ 
above a few TeV x (eV/T) is then 



_ Va' 1 



(7) 



where the flux J-'Ai— cua) is integral over the energy- 
decade of the peak region, A (eV/T) PeV < Ea < 
10A(eV/T). It is commonly assumed that this nuclear 
flux results from continuous trapping of the diffuse cos- 
mic ray flux by diffusion in a milligauss magnetic field. 

Putting into Eqs. ^ and ([7]) the parameters for the 
Cygnus 0B2 region, i.e., A^ = (56/^1) x 10^^ cm, i?sB = 
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10 pc, and d = 1.7 kpc, one obtains = 2 x 10~ A^^. 
Thus, an accumulated nuclear flux within Cygnus 0B2 
of //cm^/s in the PeV region gives a > TeV 7-ray flux at 
Earth of ~ lOA/ x Fcrab(eJ;^^ > 1 TeV). For example, 
an iron flux T^q = 6 x 10^^/cm^/s above a TeV leads to 
a 7-ray flux at Earth of about 3% of the Crab, just at 
the level measured by the HEGRA experiment for 7-ray's 
from the Cygnus 0B2 direction [l^ . In our accompany- 
ing paper [7| we present a more detailed calculation of 
the A*-process for this starburst region, and show that 
-^56 ~ 10~''/cm^/s is a credible 1% of the kinetic energy 
budget of Cygnus 0B2. 

The 7-ray energy spectrum remains to be discussed. In 
the rest frame of the excited nucleus, the photon is emit- 
ted isotropically and nearly monochromatically. There- 
fore, in the lab frame, the 7-ray-spectrum is nearly flat 
between and 2 Ta on a linear scale. The power 
spectrum and integrated power rise as eJ^;^^ and (e!j;^^)^. 



respectively, peaking near e 
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10TeV(eV/T), as explained in the paragraph below 
Eq. Q. 

In detail, the photon spectrum is obtained by 
replacing the approximate Eq. fT]) with an inte- 
gral over {dFA/dE\^^) x (Aj|^)"\ with measure 
dE^^ dcos9-y de^AB ^[gLAB_p^ ^dxn (i+cos 61-.,)], where 
9^ is the angle in the nucleus rest frame between the 
isotropically-emitted photon and the boost direction. Af- 
ter assuming a power-law (with spectral index a) for the 
nuclear spectrum, and integrating over dcosOj, one ar- 
rives at 



de^AB 
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where ri^^ is the true density of the thermal photons 
after spatial dilution, Eq is any reference energy, ~ 
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variable, and 



e^"" T/e5:™e!:^^" is the dimensionless energy-scaling 
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is the scaling function. 
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Shown in Fig. 



(9) 
© is 

Xe'F{xe) [15|. The 7-ray energy is e^^^ = 10 TeV x 

(e^=^"/MeV) (e^D^/lOMeV) (T/eV)-i. We note that the 

prefactor of F{xe) in Eq. ([8]) scales in A for fixed en- 
ergy per nucleon as ct^idr ^ (eGDR)-(a-i) ^ ^l+0.21(a-l)^ 

while the position of the peak in the x^'^ F{xe) spectrum 
at Xe ~ 0.25 scales mildly as x^ cx e™^ cx A'^ '^^. 
Well beyond the peak, the x^^ F{xe) spectrum falls as 
x~^ InXe; equivalently, the 7-ray spectrum falls as ~ 

(eLAB)-4. 

The Cosmogenic A*-Process. One well-known ap- 
plication 16] of nuclear photo-disintegration occurs in 



the calculation of the propagation of extreme-energy cos- 
mic nuclei in the cosmic microwave background (CMB) 
and cosmic infrared background (CIRB). The CMB con- 
tribution is thermal, with temperature T = 2.3 x 10~^ eV, 
and Eq. ([5]) readily yields a mfp of ^ A^^ Mpc for nuclei 
with energies near Ea ~ 4 A x 10^^ eV. The CIRB, now 
fairly well known [l^, is much smaller than the CMB. 
Its spectrum may be approximated by one or two power- 
laws, which allows an analytic integration of Eq. ([3]). The 
result is a mfp longer than that of the CMB result, but 
relevant to nuclei with energies down to 10^^ eV. 

Of interest in our work is the 7-ray flux produced when 
the photo-dissociated nuclear fragments produced on the 
CMB and CIRB de-excite. These 7-rays create chains of 
electromagnetic cascades on the CMB and CIRB, result- 
ing in a transfer of the initial energy into the so-called 
EGRET region below 100 GeV, which is bounde d by ob- 
servation to not exceed LUcas 2 x lO^'^eV/cm^ [20|. 

Fortunately, we can finesse the details of the calcula- 
tion by arguin g in analogy to work already done. Two 
recent papers 21, have calculated the Ve flux re- 
sulting from photo-disintegration of cosmic Fe, followed 
by /3-decay of the associated free neutrons. The photo- 
disintegration chain produces one /3-decay neutrino with 
energy of order 0.5 MeV in the nuclear rest frame, for 
each neutron produced. Multiplying this result by 2 to 
include photo-disintegration to protons in addition to 
neutrons correctly weights the number of steps in the 
chain. Each step produces on average one photon with 
energy ~ 3 MeV in the nuclear rest frame. Compar- 
ing, about 12 times more energy is deposited into pho- 
tons. Including the factor of 12 relating lOj to uji,, we 



find from Fig. 3 in Ref. (21| that cosmogenic photo- 



disintegration/de-excitation energy is more than three 
orders of magnitude below the EGRET bound . This 
result appears to be nearly invariant with respect to vary- 
ing the maximum energy of the Fe injection spectrum 
(with a larger i?max, the additional energy goes into cos- 
mogenic pion production). Thus, there is no constraint 
on a heavy nuclei cosmic-ray flux from the A* mecha- 
nism. 

Conclusion. In final summary, we have presented an 
alternative mechanism (A* ) for generating TeV 7-ray's in 
starburst regions. It has a unique orphan-like signature, 
with a flat spectrum below ^ 20 TeV/ (T/eV) and a quasi 
power-law above. Since the EM and PION mechanisms 
also produce unique signatures, gamma-ray astrophysics 
benefits from a one-to-one correspondence between the 
source dynamic and the observable spectrum. 
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